clear all
path(path,'../TOOLS_01');

[parametros_malha] = le_malha('Malha_problema_direto_5_nos.msh');

n_eletrodos=parametros_malha.n_eletrodos;
n_elementos=parametros_malha.n_elementos;
n_nos=parametros_malha.n_nos;
n_objetos=parametros_malha.n_objetos;
coords=parametros_malha.coords;
topologia=parametros_malha.topologia;
topologia_eletrodo=parametros_malha.topologia_eletrodo;
idx_malha=parametros_malha.idx_malha;
idx_objeto=parametros_malha.idx_objeto;




%parametros gerais
altura=0.05; % altura da malha
rho_eletrodo=1/50; % resistividade do eletrodo
corriente=1e-3; % corrente aplicada

%parametros da malha
rho=2.0;
rho_objeto=200;


rho_malha=zeros(1,n_elementos+n_eletrodos);
rho_malha(idx_malha)=rho;
rho_malha(idx_objeto)=rho_objeto;
rho_malha(n_elementos+1:end)=rho_eletrodo;




V_objeto = p_direto(corriente,n_nos,n_elementos,coords,topologia,rho_malha,n_eletrodos,topologia_eletrodo,altura ,topologia_eletrodo);

save('lixo.mat','V_objeto')
% 
% figure(1)
% for i=1:n_eletrodos
%     plot(V_objeto((i-1)*32+1:i*32),'x'),grid
%     pause
% end


